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Abstract 

The Rouse (J. Chem. Phys. 21, 1272 (1953)) and Zimm (J. Chem. Phys. 24, 269 (1956)) treat- 
ments of the dynamics of a polymer chain are shown to contain a fundamental mathematical error. 
As a result, the oft-cited mode solutions for these models are qualitatively incorrect. Compari- 
son with the Wilson-Decius-Cross treatment of vibrational modes of polyatomic molecules reveals 
qualitatively the correct form for the solutions to the Rouse model. 
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For the past half-century, the Rouse|l|| and Zimm[2| models have provided a theoretical 
framework for the description of polymer dynamics in systems ranging from dilute solu- 
tions to the melt. The objective of this note is to identify a critical mathematical defect 
in the Rouse mode solution for these models. The nature of the failing is immediately evi- 
dent from the Wilson-Decius-Cross|3|] treatment of small- displacement vibrations of isolated 
molecules. The failure, which is fundamental and arises at a mathematical rather than a 
physical level, has as a result that the characteristic modes and characteristic frequencies 
usually identified as solutions for these models are incorrect. The following presents first the 
salient features of the Wilson-Decius-Cross and the Rouse-Zimm models, shows where their 
solutions contradict each other, and demonstrates where the Rouse-Zimm treatment went 
astray. 

Wilson, Decius, and Cross [3] treat small vibrations of an iV-atom molecule. For the 
displacements r-j(t) of the atoms from their equilibrium locations, the equations of motion 
may be written 



Here m« is the mass of atom i and the Vjj are components of the matrix of second derivatives 
(a 3iV x 3N matrix) of the potential energy. Eqn. [1] represents a set of 3N coupled linear 
differential equations. The solutions to these equations are a set of 3iV eigenvectors (repre- 
senting atomic displacements) and 3N corresponding eigenvalues A«, the Aj being squares of 
the corresponding vibration frequencies. Six eigenvectors, each having Aj = 0, correspond 
to whole-molecule translations and rotations that do not alter the relative positions of the 
atoms. The remaining 3iV — 6 internal normal modes represent the molecular vibrations. 

The degenerate zero eigenvalues create technical mathematical difficulties with solving 
the corresponding eigenvector-eigenvalue problem. Wilson, et al.[3j solved the problem by 
identifying an appropriate set of internal molecular coordinates, including bond stretches, 
bends, torsions, and out-of-plane motions, and a matrix method for transforming eqn.[T]from 
its cartesian 3iV-dimensional form to a new 3N — 6-dimensional form. In the new form, the 
basis vectors are internal coordinates, while the whole-molecule translations and rotations 
occupy a 6-dimensional subspace orthogonal to the subspace of the internal coordinates. 
The Wilson-Decius-Cross FG method thus leads to a non-singular matrix whose solutions 
are the 3N — 6 internal modes. One also implicitly obtains six eigenvectors having eigen- 
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value zero and corresponding to whole-molecule translations and rotations. For the simplest 
case of a bent triatomic molecule A-B-A, the 3N solution vectors are three orthogonal free 
translations, three orthogonal free rotations, and three internal modes qualitatively and im- 
precisely described as the bending mode, in which the A-B-A bond angle oscillates, and two 
bond stretching modes, in which the two A-B bond lengths oscillate in or out of phase with 
each other. 

Comparison is now made with the Rouse and Zimm polymer models. These models 
treat a linear chain of frictional beads ("sub-molecules"), here labeled (1,2, ... ,N). The 
coordinate of bead % is r$, the vector displacement from bead % — 1 to bead %, i = 1 being 
a special case. Because their origins move, the form a set of non-inertial coordinates, 
but effects arising from the non-inertiality will be seen to be negligible. Each bead has 
mass m and drag coefficient /. The distribution of bead-to-bead distances corresponds to 
a potential of average force whose gradients are the thermally averaged forces that beads 
exert on each other. The distribution of distances is Gaussian, so the average force between 
adjoining beads is a harmonic restoring force having a spring constant here denoted k. For 
a polymer chain in a stationary solvent, the Newtonian equations of motion corresponding 
to the Rouse model are (except for special-case beads 1 and N) written 

= -f^jp- + k(v i+1 (t) + r^t) - 2r t (t)) (2) 

At low frequencies of interest here, the inertial terms forming the LHS of this equation, 
and matching fictitious forces arising from the non-inertial nature of the coordinates, are 
negligible. If the solvent is not stationary, additional forces fu^ appear, Ui being a local 
solvent velocity at the location of bead i. From symmetry, (rj) = 0, so the do represent 
displacements of the beads from their current equilibrium positions. 

Several simplifications are now applied. The bead-bead average force represents a bond 
stretch. The angles between pairs of adjoining bonds do not affect the average forces. Beads 
i and i ± n, n > 1 have no direct interactions. Therefore, the x, y, and z displacements of 
the different atoms are seemingly uncoupled, allowing eq. [2] to be reduced to a set of N — 2 
equations for the x coordinates of beads 2, 3, . . . , N — 1, 

f-^ 1 = k(x i+1 (t) + x^(t) - 2xi(t)), (3) 

plus two-special-case equations for the x coordinates of beads 1 and N, and two correspond- 
ing sets of N equations for the y and z coordinates. The Zimm model adds to the forces 



on bead % the hydrodynamic forces ^ • Ty • Fj due to mechanical forces Fj exerted on the 
solvent by the other beads j, Ty being the pre-averaged Oseen tensor in which Fj and 
Ty • Fj are parallel. 

Eq. [3] can be converted to a Langevin-type equation by adding a random thermal force 
Ti{t) on each bead. Alternatively, after adding thermal forces —ksTV^, Zimm converted 
eq. [3jto an equation for the distribution function \l/(ri, . . . , r^v) of chain ends. These changes, 
while central to applications, are not relevant here. 

Equations [3J have for eigenvalues a single zero, corresponding to a mode in which all iV 
beads are displaced the same distance in the x direction, and N — l non-zero, non-degenerate 
relaxation times corresponding to N — 1 internal modes for the Xj. The relaxation times are 

r = I (£) 

n 8ksm 2 {mr/2N) K ' 

for n G (1, 2, ... , N—l). The displacements of the atoms are determined by the amplitudes 
C n of the normal modes via 

x l = C + 2 N ^C n cos( 7in{l - 1/2) ), (5) 

71=1 ^ ' 

Cq being the center-of-mass location, and with corresponding solutions for motions in the y 
and z directions. The relative atomic displacements in a single mode n are found by taking 
C n = 1 and all other Cj, j ^ n, to be zero. 

The above coupled differential equations are separable, so to solve their spatial parts 
(the RHS of eqs. CD and [3]) it does not matter that their temporal parts (the LHS of the 
eqs. [1] and [3]) differ. Furthermore, from both a mathematical and a physical standpoint, 
the spatial parts of the Wilson-Decius- Cross and Zimm pictures, the RHS of eqs. CD and [3l 
respectively, are the same. In each case, one has a set of coupled linear differential equations 
with constant coefficients. In each case, the linear couplings are restoring forces that are 
linear in the distances between pairs of moving objects. It is of no consequence that in one 
case the moving objects are called "atoms", and in the other the moving objects are called 
"beads". It is of no consequence that in one case a coordinate is a displacement from an 
absolute equilibrium location, and in the other case a coordinate is a displacement from the 
current location of the prior bead in the chain. 

The spatial parts of the Wilson-Decius-Cross and Rouse problems are thus identical in 
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certain mathematical properties, so therefore the aspects of their solutions that are deter- 
mined by those properties must also be the same. 
Unfortunately, they are not. 

The remainder of this paper exhibits for a simple model case the contradictions between 
the Wilson-Decius- Cross and Rouse solutions, and explains why the Rouse mode solutions 
are incorrect. 

As a simple model, consider a bent triatomic molecule, such as water. A single H-O-H 
molecule has six modes (three translations and three rotations) that do not displace the 
atoms with respect to each other, so that they have Aj = 0. The H-O-H molecule also 
has three internal vibrational modes with non-zero eigenvalues, namely a symmetric stretch 
mode, an antisymmetric stretch mode, and a bond bending mode. 

Contrast the nine water modes with the modes of a three-bead Rouse model polymer, 
which may be obtained from eq. [5] by setting each of the Cj seriatim to a non-zero value. For 
the x coordinates, the polymer has one mode of zero frequency, corresponding to uniform 
translation, one mode in which the first and third beads have equal and opposite displace- 
ments parallel to the x-axis while the center bead is stationary, and one mode in which the 
first and third beads move in the same direction along to x-axis while the center bead moves 
twice the distance in the opposite direction. The modes for the y and z coordinates are the 
same (except for direction) as the x modes. The six stretch modes of the Rouse model are 
therefore trebly degenerate. The Rouse model three-bead polymer has only three modes 
having eigenvalue zero, and six modes with non-zero relaxation times. 

The disagreements between the Wilson-Decius-Cross and Rouse model normal modes 
are substantial. Does a three-mass object with purely internal forces have six modes with 
eigenvalue zero, as predicted by Wilson, Decius, and Cross, or only three, as predicted by 
the Rouse modes? Do the internal modes have three distinct eigenvalues, or do they have 
two trebly degenerate eigenvalues? The answer is supplied by classical mechanics, which 
proves that purely internal forces that act along the lines of centers between masses within a 
body cannot affect the angular momentum of that body. The correct solution to the Rouse 
model, in addition to its three translational modes having Aj = 0, must also have three 
rotational modes not seen in eqs. H] and [5] having A, = 0. 

Wherefrom comes the disagreement? Whole body rotations do not change bond lengths 
because they make simultaneous changes in bead positions along several axes. To keep bond 
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lengths constant during rotation, changes in bead locations along the a;- axis are compensated 
by simultaneous changes in bond locations along the y and z axes. The Rouse calculation 
breaks bead motions into three separate subspaces corresponding to the three coordinate 
axes, so that it has no way to generate the coupled displacements in different directions that 
represent rotations without changes of bond length. 

Mathematically, what has gone wrong is that the Rouse model solution applies stan- 
dard eigenvalue-eigenvector methods to obtain the modes and relaxation rates of a singular 
matrix, namely the singular matrix Vij of coupling coefficients. The standard eigenvalue- 
eigenvector method gives invalid answers if the matrix of coefficients is singular, precisely 
as seen here. The so-called Rouse modes, obtained by applying the standard eigenvalue- 
eigenvector method under circumstances such that the standard method is invalid, are not 
the actual modes of the Rouse model. 

It might be proposed that the number of Rouse modes is actually extremely large, so 
even if only a few of them were incorrect the consequences might not be serious. However, 
the zero-relaxation-rate modes that the Rouse modes fail to generate are the whole-body ro- 
tations. Whole body rotations are the dominant modes for viscous dissipation and dielectric 
relaxation, and cannot be neglected. 

It might be possible to apply some process analogous to the Wilson GF method to solve 
the Rouse model, at least for a chain in a particular given conformation. The potential energy 
matrix and set of restoring forces proposed by the Rouse model are purely and entirely a 
set of stretch restoring forces lying along the N — 1 bonds connecting the N beads of the 
chain. There are no forces analogous to the Wilson-Decius-Cross three-atom and four-atom 
restoring forces that resist changes in bond angles, bond torsion angles, and bond out-of- 
plane displacements. If there are no restoring forces that resist, e.g., changes in the angle 
between neighboring bonds, then modes that substantially act by changing those angles 
have eigenvalue zero. For the three-bead model system, there is no restoring force in the 
bead-bead-bead bond angle, so the Rouse mode directly analogous to bond bending should 
have Aj = 0. Should a GF-type solution route be attempted, the qualitative form of the 
solutions is immediately evident. The potential energy can be written in terms of precisely 
N — 1 distance coordinates that determine all the non-zero contributions to the potential 
energy, so such a method would decompose the RHS of eq. [3] into a 2N + 1 dimensional 
subspace having only zero eigenvectors, and an N — 1-dimensional subspace whose N — 1 
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modes are the true modes of the Rouse model. 

It is unclear that a Wilson-Decius-Cross decomposition is the most profitable way to 
tackle this problem. As an alternative, one might begin by calculating the contribution to 
polymer diffusion, viscosity, rotational diffusion, or dielectric relaxation of only six modes, 
namely the three translational modes and the three whole-body rotational modes. The 
internal modes (necessarily orthogonal to the six chain modes being included) might then 
be treated as providing perturbations. The original calculation of rj and D s by Kirkwood 
and RisemanlJ] implicitly uses this approach, while terming internal modes "fluctuations". 
In this alternative approach, one remains in Cartesian coordinates, thereby permitting a 
straightforward analysis of hydrodynamic interactions between polymer coils, as seen in 
Refs. 5 andO. 

To return to the primary conclusion, the standard mode solution to the Rouse model 
is incorrect. The standard solutions fail to generate the whole-body chain rotations that 
make the dominant contributions to viscosity and other polymer solution properties . The 
terms are lost for fundamental mathematical reasons, not because the models are physically 
inadequate, with the consequence that calculations based on the mode solution may be 
invalid. 
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